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Summary 

The first-order equations for the conservation of mass 
and momentum in differential form are combined for an 
ideal gas to yield a single second-order partial differential 
equation in one dimension and time. Small-perturbation 
analysis is applied. A Fourier transformation is per- 
formed that results in a second-order, constant-coefficient, 
nonhomogeneous equation. The driving function is taken 
to be the source of combustion noise. 

A simplified model describing the energy addition via 
the combustion process gives the required source infor- 
mation for substitution in the driving function. This 
enables the particular integral solution of the nonhomo- 
geneous equation to be found. This solution multiplied 
by the acoustic pressure efficiency predicts the acoustic 
pressure spectrum measured in turbine engine 
combustors. The prediction was compared with the 
overall sound pressure levels measured in a CF6-50 
turbofan engine combustor and found to be in excellent 
agreement. 


Introduction 

The reduction of aircraft noise through the use of high- 
bypass-ratio turbofan engines has revealed another noise 
source for concern: core engine noise. Core engine noise 
is defined as the noise produced by the gas generator- 
turbine-tailpipe combination used to provide shaft power 
to drive the fan in a turbofan engine. Reference 1 sum- 
marizes various combustion noise prediction methods for 
aircraft engines available before 1974. The methods 
include both engine and combustor rig prediction. 

Since the publication of reference 1 a data bank from 
various combustion noise experiments has been created 
by NASA (refs. 2 to 9). Additionally, a number of engine 
tests have been performed (refs. 10 to 19) to determine 
the effect of core engine noise on the far -field observer. 
The acoustic power determined from the data in refer- 
ences 2 to 9 has been correlated with engine operating 
parameters in reference 20. A number of theoretical 
investigations are reported in the literature: reference 21 
sets forth several possible combustion noise source terms 
to be evaluated; reference 22 gives a correlation of data, 
based on theory, for the acoustic power in a combustor 
with engine operating parameters; reference 23 gives 
scaling laws for combustion noise radiating from open 
flames that are based on a theory evolved by using the 
chemical reactions; reference 24 calculates the 1 /3-octave 
spectrum for a propane-air open flame; and reference 25 
points out the need to know the heat release distribution 
in the combustor before the dynamic design can be 
performed. References 21 to 24 have focused attention on 
combustion noise theory and have presented a funda- 


mental, though at times complex, picture of possible 
combustion noise mechanisms. 

In an effort to understand the dominant combustion 
noise mechanism and to provide a theoretically based 
prediction that allows insight into the noise-generating 
parameters in ducted combustors, a simple combustion 
noise model is postulated in this report. The classical 
small-perturbation approach is used to determine the 
acoustic source pressure spectrum and the parameters 
affecting the sound pressure. An acoustic pressure 
efficiency is introduced in the theoretical expression for 
the acoustic pressure to account for the inefficiency of 
converting turbulent energy to acoustic energy by the 
combustion process. The symbols used in this report are 
defined in appendix A. 


The Governing Differential Equations 
and Assumptions 

The general laws of conservation of mass and momen- 
tum and the assumptions used to limit the analysis to a non- 
viscous, adiabatic, perfect gas flowing through a constant- 
area duct are given here. The second-order partial 
differential equation governing the noise generated in a 
one-dimensional duct by the combustion process can be 
derived from the laws of conservation of mass and 
momentum. These laws given in equation form, taken 
from reference 26, are 

Continuity: 


^ _ dpV 

dt dx 


( 1 ) 


Motion: 


dpV ^_dpV^ _dP 

dt dx dx ' ^ 

For this analysis the following assumptions have been 
made: 

(1) Turbulent shearing stress is negligible (i.e., the 
fluid is inviscid). 

(2) Body forces due to gravity are negligible (i.e., 
pg=0)- 

(3) Heat transfer across boundaries is negligible (i.e., 
(VQ^O). 

(4) Changes in the flow area are negligible. 

(5) The fluid acts as a single-component gas (i.e., the 
effect of fuel mass in the mixture is negligible). 

(6) The fuel contributes energy to the air and other- 
wise has a negligible effect on fluctuating pressure. 



(7) The gas follows the perfect-gas law. 

(8) Turbulence is dominant over acoustic fluctuating 
quantities in the combustor. 

(9) Energy addition begins at the flame front. The rate 
of energy addition decreases exponentially with axial 
distance. 

(10) The Mach number in the combustor is much less 
than unity. 

(11) The combustor is located in an infinite pipe; 
hence, no acoustic reflections are present in the 
combustor or pipe. 

Source Model for Combustion Noise 

The model for combustion noise considers a one- 
dimensional duct flow (fig. 1). The flow is composed of a 
mixture of fuel and air. Turbulence is generated upstream 
of the flame front by tabulators. The flame front is 
stabilized aerodynamically at a point just downstream of 
the turbulators but does not touch the turbulators. The 
source region of combustion-related noise starts at the 
flame front. To simplify the calculations, the flame front 
is considered to be in a plane perpendicular to the flow 
direction. The regions upstream and downstream of the 
source region have constant though different properties 
and are infinite in length. 

The objective of the present work is to describe the 
pressure fluctuations at the combustor that drive the 
acoustic system upstream and downstream of the 
combustor. 



Figure 1 Combustion noise source model. 


Acoustic Source and Wave Equations 

The acoustic wave equation describing the sound prop- 
agation in a duct is derived in appendix B by using the 
conservation equations (eqs. (1) and (2)) and assumptions 
1 to 7. The resulting wave equation (eq. (Bll)) may be 
applied to the duct either upstream or downstream of the 
combustion zone since the properties in either direction 
from the source are constant though different from each 
other. The wave equation (eq. (Bll)), simplified by 
assuming constant mean properties in the region external 
to the combustor, is written here as 

1 d2p d2p _ ( p ^ ,^d2pFv , d2dV2 

^ W '^~d(2\c2 ) dx2 dx2 

In the derivation given in appendix B it has been 
assumed that the turbulence both upstream and down- 
stream of the combustor is smaller than the acoustic 
fluctuations. It has also been assumed, for purposes of 
calculating the acoustic source in the combustor, that the 
turbulence generated in the combustor is much larger 
than any acoustic waves that pass through the com- 
bustion region. 

The acoustic pressure is usually given in terms of finite- 
bandwidth, mean-square pressure fluctuations as a 
function of frequency. This is accomplished by taking the 
Fourier transform of equation (3), which removes time as 
a variable and thus simplifies the method of solution. The 
Fourier-transformed equation (appendix B, eq. (B19)) in 
terms of wave number is 

^^+k2p^=- <f!„^^{x,ui) (4) 

n = i 

where and are shown in appendix B (eqs. 

(B13b), (B14b), and (B15b)) to be the entropy, cold-flow, 
and direct-combustion noise sources. 

Solution of Nonhomogenous 
Wave Equation 

By using the methods given in reference 27, the 
complete solution of the homogeneous part (i.e., with 
^ sO) of equation (4) is given in appendix B as 

Po>,o = Cie‘>^+C2e-‘>^ (5) 

The acoustic pressure generated by the combustor 
installed in an infinite pipe can be calculated with the help 
of Green’s function (ref. 28) as given in appendix B. For 
the left-running wave the acoustic pressure from 
appendix B (eq. (B22)) is 


2 


where 

<Pa= -(<Pl, a + <P2,u + <P3,w) (7) 

and 


An additional source results from If the entropy 
fluctuations resulting from unsteady combustion produce 
fluctuating densities and these fluctuating densities pass 
through the combustion zone (i.e., the region having a 
nonzero Laplacian of the mean velocity), euiother noise, 
also termed entropy noise, is generated. For purposes of 
the present investigation this entropy noise is also 
neglected. 


"•“■■"'(a'*") 




fV 32» 


A dx2 


_JV d2 /6„ \ 


-U>2pS' 

Cp 


( 8 ) 

(9) 

( 10 ) 


A similar result is given in appendix B for the right- 
running wave. 


Physics of Noise Generation Process 

Examination of the three acoustic source terms, 

and ^3^^ given by equations (8), (9), and (10), 
respectively, shows that each is a function of frequency. 
The first acoustic source term is a function of the 
turbulent fluctuating pressure and density. For a caloric- 
ally and thermally perfect gas with sonic velocity Q, 
can be expressed in terms of the fluctuating entropy (eq. 
(8)). The conclusion is drawn then that the first source 
term is a result of unsteady burning that may be the 
result of nonuniform fuel droplet size or mixing. In this 
work the entropy fluctuations are neglected. 

The second acoustic source term (^)) is given 

by the second derivative of the fluctuating velocity with 
respect to axial distance x, (The second derivative with 
respect to x is hereinafter called the Laplacian.) 
Turbulence, generated by turbulators in the combustor, 
is used to promote the mixing of the fuel and air prior to 
burning. The product of the mass flow rate per unit 
cross-sectional area and the Laplacian of the turbulent 
velocity describes a noise that is present even if the flame 
is off. Hence, ^2,0? appears to be a cold-flow noise source. 

The third acoustic source term ^3^^^ (eq. (10)) is due to 
the Laplacian of the product of the turbulent density and 
the mean velocity. It has been assumed in this work that 
the turbulent density fluctuations in the flame zone are 
constant throughout the zone. With this assumption the 
source term becomes a function of the Laplacian of the 
mean velocity, the result of the heat energy added by the 
combustion process. This source is investigated in detail 
in order to find the scaling coefficient or correlating 
parameters applicable to combustion noise. 


Acoustic Source Term Evaluation 

Neglecting the temperature fluctuations and the 
cold-flow noise ^2,0? generated by the turbulators leaves 
the direct-combustion noise source to be evaluated. 
Equation (6) calculates the acoustic pressure due to the 
acoustic source terms. 

To evaluate (eq. (6)), the source term ^3,^0 must be 
given as a function of distance x. It is shown in appen- 
dix B that for Mach numbers much less than unity the 
velocity is directly proportional to the heat energy input 
to the air (appendix B, eq. (B36)). 

By assuming that the fuel droplet mass decreases 
exponentially with axial distance, the velocity 
distribution due to the energy released to the air can be 
determined as a function of x. The second derivative with 
respect to x is then substituted into equation (10). The 
Fourier-transformed source term from appendix B (eq. 
(B43)) is 


^ 3 , 0 ) = 


A 0 CpTr^o 


(11) 


Inserting equation (1 1) into equation (6), assuming that 
the fluctuating density is constant with distance through 
the combustion zone, and integrating yield the acoustic 
pressure in polar form as (see derivation in appendix B, 
eq. (B46)) 

VqHJq X e‘^kx-(T/2) + tan - •( - Ar/X)] ^ 
2 p A CpTr,o k ^l+(/t/X)2 ^ ^ 

where 0<jr<L^. Equation (12), with dimensions of 
pressure per unit angular frequency, is used for calcu- 
lating the overall sound pressure level and the narrow- 
band acoustic spectrum. 

In deriving equation (12) the conversion of turbulent 
energy to acoustic energy was assumed to be complete. 
This, in reality, cannot be accomplished as evidenced by 
the fact that turbulence does exist downstream of a 
combustor. To account for this, an acoustic pressure 
efficiency rip is introduced and defined as the ratio of the 
experimentily measured acoustic pressure generated by 
the combustor to the acoustic pressure that would be 
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and (16) as 


i llliiilii 


generated if all of the turbulent energy had been con- 
verted to an acoustic pressure From appendix B (eq. 
(B45a)) rip is given by 



(13) 


Thus the acoustic pressure efficiency can now be deter- 
mined for any specific combustor. The value of rip will 
probably depend on the combustor type and 
configuration. Therefore it is necessary that rip be 
determined for a number of combustors so that the 
proper value can be selected when making predictions of 
combustor pressure fluctuations. 


Overall Sound Pressure Level 

An expression for the overall sound pressure level can 
be obtained by multiplying equation (12) by the acoustic 
pressure efficiency rip, squaring, and then integrating 
over the frequency range of interest. The resulting 
equation for is 

(14) 

VPref/ \ 2pref P ^ CpTjfi) 

where the term 3^ is a function containing all of the 
frequency dependence and is defined as 


OASPL = 20 log 10 


VpVcC gg, w VqHJo 


(17) 


The acoustic pressure efficiency r)p can be determined 
from the measured overall sound pressure level and 
equation (17). The equation for experimentally deter- 
mining the value of rip is 


/ _ ^Pr^f^by/fE £_ ^ X lOOASPL/20 (18) 

The downstream temperature is used in calculating the 
sonic velocity C because the measurements in combustors 
are in the hot region. Acoustic pressure waves moving 
upstream from the combustor would require that the inlet 
temperature be known when calculating their acoustic 
pressure. 

The fluctuating pressure in the combustor of the 
CF6-50 turbofan engine has been measured (refs. 7 and 
8). The measured overall fluctuating pressure (assumed 
to be acoustic) in the combustor at 3.8 percent of design 
thrust has been used in conjunction with equation (18) to 
calculate the acoustic pressure efficiency, = 0.030, for 
the CF6-50 combustor. 

From this acoustic pressure efficiency the overall 
sound pressure level in the combustor was predicted over 
the engine thrust range by using equation (17). The re- 
sults are shown in figure 2. The theory, given by the solid 
line, agrees well with the measured fluctuating pressure, 
given by the symbols. Thus the theory can be used to 
predict, at a very minimum, the trends in the overall 
acoustic pressure. 

In light of this discussion it can be concluded that the 
theory predicts the trends in the combustor overall 


Equation (15) indicates that 3^2 is inversely pro- 
portional to the frequency. Therefore the magnitude of 3 
will be determined by the lower limit on frequency, 
providing that the upper limit is much greater than the 
lower limit. In determining the OASPL the combustion 
noise frequency limits are of the order of 50 to 2000 Hz, 
For a fuel mass decay constant X, defined as IrrfLij, and 
reasonable sonic velocities, equation (15) can be 
approximated by using the lower frequency limit by 


32 = 



, sec - 1 


(16) 


The overall sound pressure level (OASPL) equation can 
now be written for the combustor from equations (14) 



Figure 2. - Combustor pressure fluctuations as a function of CF6-50 
turbbfan engine thrust — experimental data and theory. 
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acoustic pressure with engine operating conditions and, 
for CF6-50 combustors and similar geometries, an 
acoustic pressure efficiency of 0.030 can be used in 
equation (17) to predict the combustor overall acoustic 
pressure. The acoustic pressure efficiency has also been 
determined for a YF102 engine combustor. This 
combustor, unlike the CF6-50, is a reverse-flow 
combustor and is much smaller than the CF6-50. Its 
acoustic pressure efficiency is 0.016 at 30 percent of 
design speed. Using the acoustic pressure efficiency of 
0.016 yields only a 0.2-dB error in predicting the OASPL 
at the 95-percent speed. 

Comparing the acoustic pressure efficiencies shows 
that either the CF6-50 combustor converts the turbulence 
to acoustic pressure more efficiently than the YF102 
combustor or that the turbulence levels in the combustors 
are not equal. 


Narrow-Band Sound Pressure 
Level Spectrum 

The theoretical sound pressure spectrum is given b> 
equation (12). Discarding the phase information con- 
tained in the exponential term and applying equation (13) 
to account for the efficiency of conversion of turbulence 
to acoustic pressure, equation (12) is written as 


r 2 p ALb _ 1 \ 

W f ^i+(/i6/C)2 

The right side of equation (19) represents the spectral 
shape of the acoustic pressure given by equation (12). 
The sound pressure level spectrum is given by 


SPL = 20 log 10 


2/?ref P ^Lfy 




c 


A/ 

( 20 ) 


This equation states that the acoustic pressure level is 
inversely proportional to the square of the frequency. 
The contribution to the spectrum of the terms contained 
under the radical is small for /< C. 

The predicted spectral shape is compared with the 
measured narrow-band spectra for the CF6-50 engine 
(ref. 7) in figure 3. The measured fluctuating pressure 
level in the combustor of the CF6-50 turbofan engine 
operating at 3.8 percent of design thrust is shown in 
figure 3(a). Also shown in figure 3(a) is a plot of equation 
(20) with the level matched to the measured spectrum at a 
frequency of 2000 Hz. The acoustic pressure efficiency is 
determined from equation (20) and is given by 


At 2000 Hz the measured SPL is equal to 1 10 dB and the 
narrow-band spectral acoustic pressure efficiency, com- 
puted by using equation (21), is 

t;^ = 4.5x10-2 

when is assumed to be 0.3. 

The trends in the spectral shape of both theory and 
experiment agree very well at the 3.8-percent-thrust point 
given in figure 3(a). The comparison at 99.8-percent 
thrust shown in figure 3(b), however, is not as good. The 
theory given for the 99.8-percent-thrust point uses the 
acoustic pressure efficiency obtained from the 
3.8-percent-thrust point, namely 7 ip = 4.5 x 10 ' 2 . 
Examining figure 3(b) shows that the general trend is 
correct but that large oscillations in the spectrum are not 



(a) 3. 8 -Percent thrust. 

(b) 99.8-Percent thrust. 

Figure 3. - Comparison of theory to measured sound pressure level 
spectrum — CF6-50 turbofan engine combustor. Acoustic pressure 
efficiency, 77 ^, 4.5 x 10~^; dimensionless fluctuating density, 
hjp, 0.3. 
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predicted by equation (20). Acoustic reflections from the 
turbine may account for these oscillations. 

From this discussion it is concluded that the spectrum 
predicted by equation (20) agrees substantially with the 
experimental measurements made in an operating engine. 
It can also be concluded that the simplifled combustion 
noise theory given herein substantially predicts the trends 
in the acoustic pressures generated by the combustion 
process of turbine engine combustors. However, addi- 
tional information on the turbulence in the combustor is 
required to evaluate its contribution to the combustion 
noise spectrum. 

1/3-Octave-Band Sound 
Pressure Spectrum 

The 1 /3-octave-band spectrum of the combustor noise 
can be calculated for a simplified spectrum by neglecting 
the term under the radical contained in equation (20) and 
integrating over the bandwidth: 

ASPL^lOlogjo^ (22) 

Jc 

This equation has been plotted in figure 4 with its level 
matched, as for the narrow-band spectrum, at the 
2000-Hz band for the CF6-50 operating at the 
3. 8-percent -thrust point. The 1/3 -octave-band data given 
by the symbols match the theoretical trends for fre- 
quencies above 100 Hz. At frequencies less than 100 Hz 
the data are significantly overpredicted by the theory. 

As has been discussed, the input to the theory 
considered the turbulence spectrum to be constant. The 
eddy size in the combustor must be limited to the height 
of the combustor. If the velocity is of the order of 18 
m/sec (60 ft/sec) and the annular combustor is approx- 
imately 0.15 m (0.5 ft) high, the frequency of the eddy 
passing a point must be no less than 120 Hz. This may 



Figure 4. - 1/3-Octave-band fluctuating pressure level spectrum shape 
for CF6-50 turbo fan engine combustor. Thrust, 3.8 percent of 
design. 


account for the discrepancy between the low-frequency 
theory and the data. This situation illustrates the need for 
turbulence measurements in combustors and the resulting 
spectral information. 


Concluding Remarks 

The theory derived herein can be applied, in modified 
form, to the noise generated by the turbine at low 
frequencies. This is apparent by inspection of the third 
source term. This source term indicates that any fluc- 
tuating quantity experiencing a mean flow acceleration 
(combustor) or deceleration (turbine) will produce a low- 
frequency sound (eq. (6)). For the turbine case the source 
term ^3^^^ becomes a function of the product of the 
fluctuating density and the Laplacian of the mean 
velocity through the turbine rotor blade. 

The acoustic pressure generated by the combustor can 
be expressed in several different forms. For instance, the 
acoustic pressure can be said to be directly proportional 
to the square of the volumetric heat release rate, the 
turbulence density fluctuations, and the mean inlet 
velocity and inversely proportional to the entrance 
enthalpy and frequency. Because other variations in the 
governing equations are easily derived, it is under- 
standable that the literature contains so many expressions 
for the acoustic pressure generated by the combustion 
process. The conclusions drawn herein are the result of 
the expression for the acoustic pressure generated by the 
interaction of the turbulence with the Laplacian of the 
mean energy created in the combustor. The existence of 
cold-flow noise and time-varying combustion has been 
noted, but they have not been evaluated. 


Conclusions 

By using the theory developed herein, expressions for 
the acoustic pressure generated by the combustion 
process have been derived. It has been shown that the 
overall acoustic pressure in large turbofan engine 
combustors can be predicted over the range of engine 
operating speeds, providing that the acoustic pressure 
efficiency, defined herein, can be determined. The 
following conclusions are based on the theory: 

1. The major source of combustor pressure fluc- 
tuations is the interaction of the turbulence and the mean 
internal energy (i.e., temperature) additions in the 
combustor. 

2. For a typical large turbofan engine the ratio of the 
measured overall acoustic pressure in the combustor to 
the theoretical acoustic pressure (obtained by assuming 
that all of the turbulent energy is converted to acoustic 
pressure) is a constant with a magnitude of the order of 
0.030. For smaller reverse-flow combustors the acoustic 
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pressure efficiency has been found to be of the order of 

0 . 020 . 

3, The turbulence-flame-generated noise is directly 
proportional to the square of the sonic velocity, the mass 
flow rate per unit cross-sectional area, the ratio of heat 
energy added per unit mass of air to the inlet enthalpy, 
the inlet velocity, the combustion efficiency, and the 
turbulence intensity. It is inversely proportional to the 


square of the frequency and the burning length of the 
combustor. 

4. Because combustor size limits turbulence scales, 
larger scale low-frequency turbulence is suppressed. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 1, 1983 
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Appendix A 
Symbols 


A cross-sectional flow area, m2 

C sonic velocity, m/sec; constant of integra- 

tion 

Cp specific heat at constant pressure, J/kg K 
Co flow or sonic velocity used in deriving 

generalized acoustic wave equation, m/sec 
Cj,C 2 constants used in complementary function 

for solving differential equations 
/ frequency, Hz; function 

lower limit on frequency, Hz 
upper limit on frequency, Hz 
/o ratio of mass of fuel to mass of air 

Af frequency bandwidth 

g acceleration of gravity, m/sec2 

Hj total enthalpy, equals heat energy released 

by combustion, J/kg (cal/kg) 

//y heating value of fuel, J/kg (cal/kg) 

/ 

J mechanical equivalent of heat, J/kg m 

k acoustic wave number, w/Cq, 1/m 

burning length (i.e., axial distance 
measured from flame to point in com- 
bustor where chemical reaction ends), m 
M Mach number, V/C 

SPL sound pressure level, dB (re 20 ptPa) 

OASPL overall sound pressure level, dB (re 20 /^Pa) 
P instantaneous or mean static pressure, Pa 

p time-fluctuating part of instantaneous 

pressure, Pa 

/?ref reference pressure used in calculating 

sound pressure level, 20 pPa 
Fourier -transformed fluctuating pressure. 

Pa sec 

Po),L acoustic pressure for left-running wave 

(i.e., xr<0), Pa sec 

Po 3,0 acoustic pressure outside source region 

(i.e., ^<0 and x>L^), Pa sec 
p^ acoustic pressure generated in combustor. 

Pa sec 

Q heat energy, J/kg (cal/kg) 


R universal gas constant for air, m/K 

S' fluctuating component of entropy, J/kg K 

T mean temperature, K 

T' instantaneous or fluctuating temperature, 

K 

t time, sec 

U instantaneous or mean internal energy, 

J/kg 

V mean velocity, m/sec 

V instantaneous or fluctuating velocity, 

m/sec 

W mean mass flow rate, kg/sec 

X axial distance measured from flame front, 

positive in direction of flow, m 
ASPL sound pressure level difference 

5 time-fluctuating component of density, 

kg/m3 

7}c combustion efficiency 

acoustic pressure efficiency 

6 phase angle, deg 

X fuel mass decay constant, assigned a value 

of 2ic/Lt, 

I axial distance in combustor measured from 

flame front location, m 

p instantaneous or mean density of fluid, 

kg/m 3 

Fourier-transformed source terms (eq. (4)) 
0 ) angular frequency, 27 t/, rad/sec 

Subscripts: 

c 1 /3-octave center frequency 

L left running 

n acoustic source numbers 

R right running 

T total conditions 

CO Fourier -transformed quantity; spectrum 

0 upstream axial location 

1 axial location downstream of combustion 

zone; first source term 

2 second source term 

3 third source term 
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Appendix B 

Derivation of Governing Differential Equation Neglecting Heat Transfer, 

Body Forces, and Shearing Stresses 


Under the assumption of negligible shearing stress, 
heat transfer, and body forces the time-dependent con- 
servation equations in one dimension can be written as 

Continuity: 


parameters in equation (B6) are instantaneous values of 
P, p, and F. They are assumed to be the sum of a mean 
value and a fluctuating component. In equation form 
they are defined as 


dp _ dpV 
dt dx 

Motion: 


(Bl) 


P=P(x)+p(x,0 
p=p(x)-\-8{xJ) 
V= V{x)^v{xj) 


(B7) 

(B8) 

(B9) 


d^ _ d^ 
dt dx dx 

Taking the d/dt of equation (Bl) gives 


Substituting equations (B7), (B8), and (B9) into equa- 
tion (B6), dropping the terms containing time derivatives 
of mean quantities, and expanding the products yield an 
equation in both mean and fluctuating quantities. 
Equating the zero-order terms (mean quantities) yields 


52 p ^ d'^pV 


dip a2pF2 ^ 

(B3) 'W- “ax2“ - 


(BIO) 


Taking d/dx of equation (B2) gives 

dipV _ d2pV2 d2p 

dxdt dx^ dx^ 

Adding equations (B3) and (B4) gives 

d2p _ a2pj/2 ^ a2p ^ 
dt2 dx2 dx2 


Equating the first-order terms yields 


1 d^p d^p _ ^2 

(B4) q 


d2p d2p d2 / p \ 


For steady flow, continuity gives a relationship between 
the mean values of p and V as 


(B5) 


pV=- 


(B12) 


Rearranging terms and adding 

to both sides of equation (B5) give 
1 32p_^_ ^ 1 32/> S^pV2. 

q aa dx2 ~ q dfi 3x2 


(B6) 


The fluctuating quantities are assumed to be the result of 
the turbulence generated by the turbulators just upstream 
of the flame front. 

Ribner (ref. 29) in discussing equation (Bll) follows 
Lighthill’s interpretation of the source terms on the right 
side of equation (Bll). Lighthill (ref. 30), neglecting the 
flrst term on the right, uses the turbulent stresses to 
formulate his jet noise theory. The source terms on the 
right side of equation (Bll) are, for simplicity, defined as 
follows: 


where Cq can be defined generally, at this point in the 
derivation, as the velocity of propagation of either a 
sonic wave or a disturbance such as turbulence. The 



(B13a) 
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or for Co deflned as the acoustic velocity 



as the fluctuating entropy source 


d^pVv 

or by using equation (B 12 ) as 

W d^v 
A dx2 

_d25V2 
~ dx2 


The first source term <pi can be written in terms of the 
fluctuating entropy by assuming that the gas is both 
calorically and thermally perfect and that the pertur- 
(B13b) bations in pressme p and density 6 are propagated at the 
speed of sound. The result of these assumptions is 
equation (B13b). Although the mean density and specific 
heat are not functions of time, it is likely that uneven 
combustion will occur and give a fluctuating heat release 
or, equivalently, a fluctuating entropy. This could be due 
(B14a) to nonuniformity of fuel droplet size or the time-varying 
mixing process in the combustor. The first source term 
thus accounts for the noise source associated with the 
fluctuating burning process as controlled by the mixing 
process in the combustor. Although the fluctuating heat 
release is recognized as a potential noise source, it is not 
(B14b) thought to be very important (ref. 5). 

Entropy noise also results when a fluctuating quantity 
passes through a gradient in a mean quantity and may, if 
(B15a) desired, be included in v ’3 by using the perfect-gas 
relationship 


or again by using equation (B12) as 


_W 32 



(B15b) 


The last term v’s indicates that the second derivative of 
the product of the fluctuating density and mean velocity 
with respect to eixial distance represents a source in the 
combustor. For simplicity, the second derivative with 
respect to distance is hereinafter called the Laplacian of 
the quantity. This source term acts as a forcing function 
for the system described by the left side of equation 
(Bll). It will be shown that the particular integral 
solution of equation (Bll) with source ^3 as its forcing 
function yields an expression for the fluctuating pressure 
that agrees with measured values in a combustor and thus 
indicates that the source is dominant over the sources 
<Pi and <P 2 ‘ 

The second term <P2 considered to be the cold-flow 
noise source term. This noise is created by the turbulators 
generating a fluctuating velocity as the flow passes 
through them. The tiubulence thus produced causes 
mixing of the fuel with the air as the mixture moves 
toward the flame front. Since the product of the mean 
density and velocity is constant and equal to the mass 
flow rate per unit cross-sectional area, the source term 
can be written as in equation (B14), leaving the Laplacian 
of the fluctuating velocity as the measure of the source 
strength. The relative magnitude of this source term can 
be determined from measurements of the fluctuating 
pressures in the combustor during flameout or cold-flow 
tests. 


6 _p r 


p P T 


giving 


w 32 \ 


^3” A 3x2 [ 

\p t)^\ 

Equation (Bll) can now be written as 


(B16) 


(B17) 


C2 3t2 


32p 

3x2 


«= 1 


(B18) 


Taking the Fourier transform of equation (B18) converts 
the equation from the time domain to the frequency 
domain. In terms of frequency then 


3x2 


3 


/!= 1 


(B19) 


where k=(j>/Co. The solution to the homogeneous part of 
equation (B19), that is. 


3x2 


+k2p^s0 


(B20) 


is 
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(B25) 


(B21) 

Since the combustor is assumed to be located in an 
infinite pipe, no reflected waves are present and the 
solution for the acoustic pressure in the infinite pipe 
radiated from the source region can be found by using 
Green’s function (ref. 28). For the left-running wave 


dx dx 


From steady-state motion 


dx dx 




pikx ?Liy 

2»S. 


x<0 


(B22) 


From the perfect-gas law 


(B26) 


where the integration is performed over the source region _ P 
(burning length). For the right-ruiming wave ^ ~ RT 


(B27) 


e—ikx fLi, 


x>Lt 


(B23) 


Substituting equation (B27) into equation (B25) gives 


where 


RT dx ^ R dx\Tj 


= 0 


(B28) 




(B24) 


or expanding the second term gives 


In this report, entropy and cold-flow noise and 
are neglected. Furthermore, the temperature 
fluctuations given in equation (B17) are considered as 
entropy fluctuations and are negligible as compared with 
the pressure fluctuations. These assumptions leave the 
direct-combustion noise source term ^3 in the acoustic 
pressure equations (eqs. (B22) and (B23)). To perform 
the integration required in equations (B22) and (B23), the 
dependence of the term d^/dx^ {d/p V) on axial position 
must be known or assumed. As discussed later, 6 /p is 
assumed to be constant. To get the velocity-distance 
relationship, we first obtain the relationship between the 
velocity and internal energy of the gas at low Mach 
numbers. It is assumed that the chemical reaction rate is 
much greater than the rate of mixing of the fuel with the 
air in the combustor. This assumption appears to be well 
justified because chemical reaction times for well-mixed 
reactants are of the order of microseconds and the mixing 
of fuel with the surrounding air appears to be of the order 
of milliseconds. Hence, the combustion process is rate 
controlled by the mixing process (ref. 31). The mixing 
process is represented by an exponential decay of the fuel 
droplet diameter with axial distance Jt. This assumption 
permits the calculation of the required velocity 
distribution with distance for use in equations (B22) and 
(B23). 

The derivation using the mean-valued parameters is as 
follows: 

From steady-state continuity 


P dV 

i 

^ / 


.pdT 

RT dx 

Rn \ 



or 





(dP 

P dT^ 

\ —0 


<dx 

T dx, 

|-0 


(B29) 


(B30) 


Next, substituting equation (B27) into equation (B26) and 
rearranging give 


dP ^ PVdV 
dx^ RT dx 


(B31) 


Substituting equation (B31) into equation (B30) and 
rearranging yield 


dx ^ ^dx) 


(B32) 


or 


/ V2\dV_VdT 
V RTjdx T dx 


(l-A/2) 


1 ^ 
V dx 


1 ^ 
T dx 


(B33) 


(B34) 
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ill 


Because the Mach number in a combustor is usually of 
the order of 0.05, equation (B34) can be simplified as 
follows under the assumption that the Mach number is 
much less than unity: 


1 1 
V dx T dx 


or 


V=Vo 


Tt 


(B35) 


(B36) 


centerline, which indicate that the bulk of the reaction 
takes place near the flame front. Reference 31 also shows 
that the turbulence decays over a large distance as 
compared with the maximum heat release region at the 
flame front. Thus, it is reasonable to expect that the 
fluctuating density can be approximated by assuming 
that it is constant near the flame front. Of course, this is 
desirable because it simplifies the solution of the 
problem. 

This reasoning then, to some extent, justifies the 
assumption of constant fluctuating density with axial 
distance near the flame front. This allows the source ^3 to 
be written from equation (B15b) as 


At this point it is necessary to know the eixial temperature 
distribution in the combustor. The heat energy released 
by the fuel droplet with x then is given by 

Q{x) =Ht{x) =Cp[Tt(x) - Tt.q] (B37) 


_ W b(PV 
A p c&c2 


(B41) 


Taking the second derivative with respect to x of the 
expression for V(x) (eq. (B40)) gives 


and if we assume an exponential decay of the fuel droplet 
diameter, 

Ht(x) =VcHM^-e->oc) (B38) 

The temperature ratio is given by assuming constant 
specific heats as 

(B39) 

f T,0 L T^O J 

where for convenience of computation, X is defined as the 
fuel decay constant and is set equal to 27 t/L^. This yields 
a vanishing exponential at x = 

Substituting equation (B39) into equation (B36) gives 
the required velocity distribution as 

= + (B40) 


cPV_ Vor,cHjQ-K2 
dx^ CpTfo 


(B42) 


Substituting this into the equation for ^3 (eq. (B41)) gives 
the direct -combustion noise source term as 


^3,0) = 


IV 

A p CpTjQ 


(B43) 


Assuming that the mass flow rate per unit area is con- 
stant and noting that the nondimensionalized density is 
assumed to be constant in the source region, the acoustic 
pressure can be calculated from equation (B22) as follows 
for the left-running wave: 


Wb^ Kor?^/oX2 
lik A p CpTjfi 


'i^b 

e-(\+ik)^di 

0 


(B44) 


To this point in the derivation little has been discussed 
about the nondimensional fluctuating density term 6/p in 
the acoustic source term ^3 (eq. (B15)). Indeed, little is 
known about this term. However, some arguments can be 
made concerning its magnitude throughout the com- 
bustion zone. First, it is assumed that the fluctuating 
density is caused by the turbulators just upstream of the 
flame front. Next, because of the exponential heat 
addition assumption, it is possible to state that most of 
the effect of the Laplacian of the mean velocity takes 
place over a short distance downstream of the flame 
front. This assumption is further justified in reference 31 
by gas samples taken along the combustor axial 


Performing the integration and substituting the upper 
limit at ^=L\j yield zero since is approximately 

zero. Substituting the lower limit ? = 0 yields 


K ^ VqHJq X2eikx 
2 p A CpTp^o ik{\ + ik) 


(B45) 


where r\p is defined as the acoustic pressure efficiency for 
the conversion of turbulence to acoustic pressure. 



(B45a) 
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The term i]p has been incorporated into equation (B45) to 
account for the fact that only a fraction of the turbulent 
density fluctuation is converted to acoustic pressure. 
Equation (B45) converted to polar form becomes 


- VpVc K ^ ypffvfo ^ 

2 p A CpTr,o k ^l+ik/\)2 


(B46) 


The sound pressure level is 


SPL = 201og,o^^ 

2 p 


IV VqHJq \ A/~ 

^ CpTr,Q k ^i+(/t/\)2 


-20 logioJ?ref (B47) 


The rate of decay of fuel droplets as they pass through 
the combustion zone is, by assumption, proportional to 
the exponential decay factor X, defined by equation 
(B38). The value of X must be such that the exponential 
function 

f{x) 


is approximately zero so that equation (B38) will yield 
negligible heat addition at the end of the combustion 
zone, x=Liy, If the value of the exponent is 2 t, thfe f{x) 
becomes 0.0019, and this indicates that practically no 
heat energy is released past the end of the combustion 
zone. For X = 2x/L^, the frequency parameter k/\ equals 
This can be used to simplify the SPL calculation 
(eq. (B47)). For source ^ 3 ,^^ it can be shown by solving 
equation (B23) for the acoustic pressure in the right- 
running wave that the sound pressure level of the right- 
running wave equals that of the left-running wave given 
by equation (B47). However, it should be noted that the 
phase angle for the left-running wave is given by 

and the right-running wave phase angle is 
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